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Abstract 

We analytically derive a relativistic variable Eddington factor in the relativistic radiative flow, and 
found that the Eddington factor depends on the velocity gradient as well as the flow velocity. When the 
gaseous flow is accelerated and there is a velocity gradient, there also exists a density gradient. As a 
result, an unobstructed viewing range by a comoving observer, where the optical depth measured from the 
comoving observer is unity, is not a sphere, but becomes an oval shape elongated in the direction of the 
flow; we call it a one-tau photo-oval. For the comoving observer, an inner wall of the photo-oval generally 
emits at a non-uniform intensity, and has a relative velocity. Thus, the comoving radiation fields observed 
by the comoving observer becomes anisotropic^ and the Eddington factor must deviate from the value 
for the isotropic radiation fields. In the case of a plane-parallel vertical flow, we examine the photo-oval 
and obtain the Eddington factor. In the sufficiently optically thick linear regime, the Eddington factor is 
analytically expressed as /(r,/3, ^) = |(1 -h where r is the optical depth and [3 (= v/c) is the flow 

speed normalized by the speed of light. We also examine the linear and semi-linear regimes, and found 
that the Eddington factor generally depends both on the velocity and its gradient. 

Key words: accretion, accretion disks — astrophysical jets — gamma-ray bursts — radiative transfer 
— relativity 



1. Eddington Approximation 



In the moment formalism of radiation transfer in a 
static atmosphere (Chandrasekhar 1960; Mihalas 1970; 
Rybicki, Lightman 1979; Mihalas, Mihalas 1984, Shu 
1991; Peraiah 2002; Castor 2004), in order to close mo- 
ment equations truncated at the finite order, we usually 
adopt the Eddington approximation in the inertial frame 
(laboratory frame) as a closure relation. 



5'^ 
3 ' 



(1) 



where E and P*-^ are the radiation energy density and the 
radiation stress in the inertial frame, respectively. 

Similary, in the moment formalim of radiation transfer 
in a relativistic flow (Thorne 1981; Thorne et al. 1981; 
Flammang 1982, 1984; Nobih et al. 1991, 1993; Park 2001, 
2006; Takahashi 2007), we often assume the Eddington 
approximation in the comoving frame (fluid frame). 



z 



1 
1 
t 

\ 


1 
1 

i 
1 


\ 
\ 

V. 


1 

* 

,/ 










/ 

1 


t 
1 
1 


1 

% 


f 
/ 

t 


\. 




/" 

t 


'\ 

1 


1 


t 



(2) 



where Eq and Pq are the radiation energy density and 
the radiation stress in the comoving frame, respectively. 
Then, this closure relation in the comoving frame is trans- 
formed to the relation in the inertial frame (Hsieh, Spiegel 
1976; Fukue et al. 1985; Kato et al. 1998, 2008). 

As is well known, the Eddington approximation is valid 
when the radiation fields are almost isotropic. In the 
relativistically accelerating flow, however, even for a co- 
moving observer, the comoving radiation fields become 



Fig. 1. Schematic picture of a relativistic radiative flow in 
the vertical direction. The flow is accelerated in the vertical 
(z) direction, and has a velocity gradient. The dashed curves 
are one-tau photo-ovals observed by a comoving observer. 
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anisotropic. 

Let us suppose a relativistic radiative flow, which is 
accelerated in the vertical (z) direction, and a comoving 
observer, who moves upward with the flow, at an optical 
depth r (figure 1). 

At the sufficiently deep inside the flow (a dashed circle 
in figure 1), where the optical depth is very large, the 
mean free path of photons is very short. In such a case, 
within the range of the mean free path, for the comoving 
observer the flow is seen to be almost uniform; the velocity 
gradient and resultant density gradient can be negligible. 
Hence, the mean free path is the same in all directions, or 
the shape of the one-tau range, to which the optical depth 
measured from the comoving observer is unity, is almost 
sphere; we call it a one-tau photo-sphere. As a result, the 
comoving radiation fields is nearly isotropic, and a usual 
Eddington approximation in the comoving frame is valid. 

In the region where the velocity is large and the density 
is low, or where the velocity gradient is large in spite of 
the high density (a dashed oval in figure 1), within the 
range of the mean free path, the velocity gradient cannot 
be neglected, and the density is no longer uniform even 
in the comoving frame of the gas. Hence, the mean free 
path becomes longer in the downstream direction than in 
the upstream and other directions, and the shape of the 
one-tau range elongates in the downstream direction; we 
call it a one-tau photo-oval. As a result, the comoving ra- 
diation fields becomes anisotropic^ and we should modify 
a usual Eddington approximation; the Eddington factor 
may depend on the optical depth r, the flow velocity ^, 
and the velocity gradient d(3/dr in the relativistic regime. 

In addition, when the optical depth is sufficiently small 
and/or the velocity gradient is suffiently large (a dashed 
hemi-circle in figure 1), the mean free path of photons in 
the downstream direction become less than unity. Hence, 
the shape of the one-tau range is open in the downstream 
direction; we call it a one-tau photo-vessel. 

In order to obtain an appropriate form of the relativis- 
tic variable Eddington factor f{r,/3^d/3/dr)^ we thus care- 
fully treat and examine the radiation fields in the comov- 
ing frame. In this paper, at the first step, we examine 
the one-tau photo-oval and derive the relativistic variable 
Eddington factor in the linear and subrelativistic regimes. 

In the next section we examine the shape of one-tau 
photo ovals seen by the comoving observer under the 
several approximations, including linear and semi-linear 
regimes. In section 3 we analytically and numerically cal- 
culate the comoving radiation fields within the photo-oval, 
and derive the relativistic variable Eddington factor un- 
der the several approximations. In section 4, we briefly 
discuss the importance of the velocity gradient on the rel- 
ativistic variable Eddington factor. The final section is 
devoted to concluding remarks. 

2. One-Tau Photo-Oval Walls 

In this section we first consider the shape of the one-tau 
photo-oval, to which the optical depth measured from the 
comoving observer is unity, in the vertical one-dimensional 




Fig. 2. One-tau photo-oval around a comoving observer in 
the vertical (z) one-dimensional radiative flow. The comoving 
observer is located at z = zq, where the flow speed is (3 = 
In the s-direction, where the angle measured from the down- 
stream direction is ^, the mean free path is set to be Iq. 

radiative flow under the linear approximation. We shall 
also obtain the breakup condition, where a photo-oval is 
open toward the downstream direction to become a photo- 
vessel (figure 1). 

The situation is schematically illustrated in figure 2. 
In the vertical {z) flow, we suppose a comoving observer 
located at z = zq^ where the flow speed \^ [5 = {(3 is the 
speed normalized by the speed of light). In the s-direction, 
where the angle measured from the downstream direction 
is 0, the mean free path of photons is Iq. The relation 
among these quantities is 



z — Zo = s COS 6 



(3) 



The continuity equation for the present one-dimensional 
relativistic flow is 



pcjl3 = J (= const.). 



(4) 



where p is the proper gas density, 7 (= — (3'^) the 

Lorentz factor, and J the mass-flow rate per unit area. In 
addition, the optical depth r is defined by 



dr : 



-f<ipdz^ 



(5) 



where is the opacity, which we assume constant in the 
present analysis. 

In this paper we use a linear approximation for the flow 
field; that is to say, around the comoving observer the flow 
speed is expanded as 



f3-- 



dz 



{z-Zq), 



(6) 



and we assume that the velocity gradient dl3/dz\Q is con- 
stant in the one-tau photo oval range. 
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Under these situations and assumptions, we examine 
three regimes below. 

2.1. Extremely Linear Regime 

In the linear regime, the density is also expanded as 

dp 



■ Po - 



Po- 



dz 
dp d/3 
dp~d^ 
pj^dp_ 
~Y~d^ 



{z - Zo) 



{z-zo) 



(z-zo), 



(7) 



where we use the continiuty equation (4), and we also 
assume that the density gradient dp/dz\o is constant. 

Then the optical depth Ts along the 5-direction is easily 
calculated as 



Jo 
Jo 



npds 



po 



Kpolo + 1^ 



dp dp 

djS dz 



dp dp 
dp dz 
dp dp 



{z - Zo) 



ds 



scosO 



dp dz 





cosO 



ds 



Hence, the length Iq of the one-tau range (r^ 
termined by the quadratic equation. 



dp dp 
dp~d^ 



-Iq^ i^polo-l = 0. 



(8) 

1) is de- 
(9) 



As is easily seen from equation (9), if there is no velocity 
gradient, tvpolo = 1^ and the shape of the one-tau range is 
a sphere (one-tau photo-sphere). 

If the deviation from the sphere is sufficiently small 
{kpqIo ^ 1)7 the shape is approximately calculated as a 
correction factor, 



hipolo = 1 — 



dp dp 



cosO 



dp dz 
dp dp 
dp~d^ 
^ I dp dp 
npQ dp dz 

Po dp dr 



cosO 1 



cosO 
cos^ 



= 1-2!^ 
p dT 



■ 
cos^ 



(10) 



Thus, a small correction is proportional to cos6>, and one- 
tau oval is elongated toward the downstream direction. 

2.2. Linear Regime 

In general the quadratic equation (9) can be easily 
solved. The quadratic equation (9) is rearranged as 




Fig. 3. Shapes of one-tau photo-ovals in the hnear regime 
for several values of parameter a. The values of parameter a 
are 0.02, 0.1, 0.2, 0.4, and 0.5 from spherical to elongated. 
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and the required solution is 
1- VI -2a cos l9 



i^Pok ■ 
where 



(11) 



(12) 



1 dp dp 
npQ dp dz 
1 dp dp 
Po dp dr 
j^dp 
' ~p"d^ 



(13) 



In figure 3 the solutions (12) of the quadratic equation 
(11) are shown for several values of parameter a. The val- 
ues of parameter a are 0.02, 0.1, 0.2, 0.4, and 0.5 from 
spherical to elongated. For a > 0.5, in the downstream 
direction the solution disappears and the photo-oval be- 
comes open to be a photo-vessel. 

From the solution (12), the breakup condition, where 
the photo-oval is open in the downstream direction, is 
given by 1 — 2acos0 < 0, or 



dp 
~d^ 



cosO > 



JL 
272 



(14) 



This breakup condition is shown in figure 4. Below the 
curve photo-ovals are closed, whereas above the curve they 
are open in the downstream direction and the present ap- 
proximation is violated. 

Using the continuity equation (4) and the optical depth 
(5), the breakup condition (14) is also written as 



dp 
dz 



o '^J 1 
COS 6 > — ^ 

2c 







(15) 



Hence, the condition is somewhat different appearance, 
whether we consider the velocity gradient in the linear 
length z or the optical depth r. 
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Fig. 4. Breakup condition for photo-ovals in the hnear 
regime. Below the curve photo-ovals are closed, whereas 
above the curve they are open in the downstream direction. 

In any case, at this stage we demonstrate that both the 
shape of the one-tau photo oval and its breakup condition 
depend on the flow velocity and the velocity gradient. 

2.3. Semi-Linear Regime 

When the flow speed is high and the density decreases, 
the mean free path of photons becomes large, and the 
one-tau range would expand. Then, even if the flow field 
around the comoving observer is assumed to be linear in 
a manner of equation (6), there is no guarantee for the 
linearity of the density distribution, and we should de- 
termine the density by continuity equation (4). We here 
examine such a semi-linear regime. 

Using continuity equation (4), the optical depth 
along the s-direction is expressed as 

= tvpds 
Jo 

^^1 1'' ^ 



ds. 



/o P 

From equations (3) and (6), we have 



(16) 



dj3-- 



dz 



df3 
dz 



cos 6 ds, 



and the integral is transformed as 



dp 
dz 



cos 6 ^ 1^0 



/3 



df3, 



(17) 
(18) 

(19) 



where we use the assumption that the velocity gradient is 
constant. This equation (19) is analytically integrated to 
give 



1 



dp 
dz 



cos 6 - 




Fig. 5. Shapes of one-tau photo-ovals in the semi-linear 
regime for several values of f3o and —df3/dT\o. The values 
of parameters are (a) Po = 0.3 and —df3/dT\o = 0.01, 0.05, 
0.10, 0.15, 0.20, 0.25, and 0.30 from spherical to elongated, 
and (b) -df3/dr\o = 0.1 and f3o = 0.1, 0.2, 0.3, and 0.4 from 
elongated to spherical in the leftside, while 0.5, 0.6, 0.7, and 
0.9 from spherical to elongated in the rightside. The dashed 
curves are photo-ovals in the linear regime shown in figure 3. 



log 



P l + yi^ 

^oi+yr^ 



dp 

dr 



-loPo 



log 



p 1 + yr^ 
^oi+yr^ 



(20) 



Thus, the length Iq of the one-tau range (ts = 1) is 
finally determined, with the help of equation (18) and 
the definition of the optical depth (5), by the following 
equations: 



0: 



dp 



cosO - 



Po 



log 



p 1 + yr^ 
^oi+yr^ 



cos^ k^pqIq. 



(21) 
(22) 
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Fig. 6. Breakup condition for photo-ovals in the semi-hn- 
ear regime. Below the solid curve photo-ovals are closed, 
whereas above the curve they are open in the downstream 
direction. The dashed curve is the breakup condition for the 
linear regime shown in figure 4. 

In figure 5 the numerical solutions of equations (21) 
and (22) are shown for several values of Po and —d(3/dr\o; 
i.e., the one-tau length k^pqIq is plotted as a function of 
0. In figure 5a the values of parameters are = 0.3 and 
-dp/dr\o = 0.01, 0.05, 0.10, 0.15, 0.20, 0.25, and 0.30 
from spherical to elongated. When the velocity gradient 
is small, a photo-oval is spherical. As the velocity gra- 
dient increases, it elongates in the downstream direction, 
and breaks up in the downstream direction at sufficiently 
large gradients. In addition, the shape of photo-ovals are 
somewhat fat, compared with those in the linear regime 
shown by dashed curves. 

In figure 5b, on the other hand, the values of parameters 
are -dp/dr\o = 0.1 and po = 0.1, 0.2, 0.3, and 0.4 from 
elongated to spherical in the leftside, while 0.5, 0.6, 0.7, 
and 0.9 from spherical to elongated in the rightside. As 
is seen in figure 5b, a photo-oval elongates in the case of 
smaller Po. This is because, for the fixed —dP/dr\o^ the 
velocity difference from the comoving observer becomes 
relatively large, if the velocity is small. 

We can further obtain the breakup condition, where the 
photo-oval is open in the downstream direction. Under 
the present approximation, equation (21) has no solution, 
when p = 1. Or, the breakup condition is 



dp 



cosO > Po 



1 



log 



l + yl^ 

Po 



1 



(23) 



This breakup condition is shown in figure 6. Below the 
solid curve photo-ovals are closed, whereas above the 
curve they are open in the downstream direction. 

3. Comoving Radiation Fields 

We now examine the radiation fields in the comoving 
frame; the radiation emitted from the one-tau photo oval 
walls forms the comoving radiation fields at the position 
of the comoving observer (figure 7). 




Fig. 7. Radiation fields in the one-tau photo-oval around 
a comoving observer in the vertical one-dimensional radiative 
fiow. The radiation fields become anisotropic, since in general 
the emitted intensity is not uniform and it is redshifted due 
to the velocity difference. 

In the static and optically thick atmosphere, the radia- 
tion fields are isotropic and uniform. In the present mov- 
ing atmosphere, on the other hand, there are two factors, 
non- uniform intensity and redshift, which make radiation 
fields anisotropic. 

First, the radiative intensity /co in the comoving frame 
emitted from the one-tau photo-oval walls is not generally 
constant, but is a function of z (or r) and fi (= cos^), 
due to the gradient of physical quantities. If the op- 
tical depth is sufficiently large, the /i-dependence could 
be safely ignored, and the emitted intensity would be as- 
sumed to be isotropic. In general, however, the radiative 
intensity increases with the optical depth, and therefore, 
the intensity from the upstream direction is slightly larger 
than that from the downstream direction in the present 
one-dimensional vertical flow. Hence, the effect of non- 
uniformity of intensity generally acts as a force to accel- 
erate a comoving observer. 

In this paper we use a linear approximation; that is, 
around the comoving observer both the flow speed and 
intensity in the comoving frame are expanded as 



dz 

r r dl 

/ = /o + :r 

dz 



{z - Zo) 



- Zo) 



dp 
dl 



i^polo, (24) 



cosO npoloi (25) 



where we use equations (3) and (5), and we assume that 
the intensity gradient dI/dz\o is constant. 

Second, the intensity observed by the comoving ob- 
server is redshifted (Doppler shifted) due to velocity dif- 
ference between the comoving observer and the one-tau 
photo-oval walls. In the accelerating flow, where the flow 
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speed increases toward the downstream direction, the rel- 
ative velocity is generally positive (figure 7) except for 
some special direction of ^ = 7r/2. Hence, the Doppler 
shift of intensity also causes anisotropy of the radiation 
fields at the comoving observer. 

The relative speed A (3 between the comoving observer 
and the one-tau photo-oval walls is given by the relativistic 
summation law, but in the linear regime we approximate 
the usual form as 



A/3: 



(3-(3o 



cos^ k^PqIq. 



Using this relative speed, the redshift z is expressed as 
1 + A/?cos6' 



(26) 



l + z = 



(27) 



Thus, the observed intensity Ico of the comoving observer 
becomes 



dJ 



1 



( d§_ 



cosO 



(^1 - ^ hipolo^ 



(28) 



3.1. Extremely Linear Regime 

Again, we first consider the sufficiently linear regime, 
where the velocity gradient is sufficiently small. In this 
case, the observed intensity (28) is linearly expanded as 



di 



cosO tipolo 



1 



Iq dr 



cosO - 



.4^ 
dr 



.4^ 
dr 



cos k^pqIq 



(29) 



where we use equation (10) for tzpolo- 

In this extremely linear regime, the radiation energy 
density Eco, the radiative flux Fco, and the radiation pres- 
sure Pco measured by the comoving observer is analyti- 
cally calculated as 



cEc, 



CO — / -i-CO^'^ ^co 



Icod^i 
27rJ, 
47rJo ( 1 



dl 

Iq dr 
4dp 
3dT 



cos6> + 4 



dp 
~d^ 



^ 1 / 12 df3 



(32) 



As already stated, if the emitted intensity from the one- 
tau photo-oval walls is not uniform, and has a gradient to- 
ward the downstream direction, as shown in equation (29), 
this intensity gradient acts as a force to push the comov- 
ing observer toward the downstream direction, as shown in 
equation (31). It should be noted that the radiative force 
received by the comoving observer is —nFco/c. However, 
this intensity gradient does not affect the radiation energy 
density and the radiation pressure. 

On the other hand, in this hear regime the relative ve- 
locity and resultant redshift is always positive in any di- 
rection except for 6> = 7r/2, and therefore, the observed 
intensity is slightly reduced in all direction, as shown in 
equation (29). This velocity gradient affects the radia- 
tion energy density and the radiation pressure, as shown 
in equations (30) and (32), but the radiative flux. That 
is, due to the velocity gradient, both the radiation en- 
ergy density and the radiation pressure slightly decrease, 
compared with the non-relativistic case. 

As a result, in this linear regime the Eddinton factor in 
the comoving frame becomes 



Eno 3 



^ 15^ 



(33) 



Hence, in the sufficiently linear regime, where the veloc- 
ity gradient is sufficiently small, the relativistic variable 
Eddington factor depends only on the velocity gradient to 
the flow optical depth. 

It should be noted that equation (33) is also expressed 



1 



16 dp 



15 tipdz 



1 



16 c-iP dp 
15 nJ dz 



(34) 



cos^O]smede 



(30) 



Fco = / IcoCosOdQc 



1 



dl 



' 3 V J^dr 
cPco = / Ico cos^ OdQc 



(31) 



Hence, in this sense the relativistic variable Eddington 
factor depend both on the velocity and its gradient to the 
real length. 

As is seen in equations (33) and (34), the Eddington 
factor slightly decreases for the subrelativsitic flow with 
the velocity gradient. 

3.2. Linear Regime 

In the linear regime, using the solutions (12), we can nu- 
merically integrate the comoving intensity (29) to obtain 
the comoving radiation fields. Since the intensity gradient 
does not affect the Eddington factor in the present linear 
treatment, we shall drop the intensity gradient in what 
follows. Hence, the parameters are the velocity Po and its 
gradient —dp/dr\o of the comoving observer. In addition, 
the integration in the polar direction is restricted in the 
range of 6i = cos~^(l/2a) < 6> < tt, although we do not 
consider a photo-vessel in this paper. 

The numerical result is shown in figure 8. In figure 8, 
3/, the three times of the Eddington factor f{P^—dP/dr), 
is plotted by a contour map as functions of the velocity Po 
(abscissa) and the velocity gradient —dp/dr\o (ordinate). 
Reflecting the breakup condition in the linear regime (see 
figure 4), the value of 3/ quickly drops at the outer edge 
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Fig. 8. Relativistic variable Eddington factor f{(3,—d(3/dT) 
in the linear regime. The value of 3/ is plotted by a contour 
map as a function of the velocity (3q (abscissa) and the ve- 
locity gradient —d(3/dT\{) (ordinate). Reflecting the breakup 
condition in the linear regime, the value of 3/ quickly drops 
at the outer edge of the one-tau photo-oval. 
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Fig. 9. Relativistic variable Eddington factor f((3,—d(3/dr) 
in the semi-linear regime. The value of 3/ is plotted by a con- 
tour map as functions of the velocity f3o (abscissa) and the ve- 
locity gradient —df3/dT\o (ordinate). Reflecting the breakup 
condition in the linear regime, the value of 3/ quickly drops 
at the outer edge of the one-tau photo-oval. 

of the one-tau photo-oval. Within the one-tau photo-oval 
region, the relativistic variable Eddington factor is well 
fitted by equation (33), except for the boundary region. 

3.3. Semi- Linear Regime 

Similary, in the semi-linear regime, using equations (21) 
and (22), we can numerically integrate the comoving in- 
tensity (29) to obtain the comoving radiation fields. The 
parameters are the velocity and its gradient —d/3/dr\o 
of the comoving observer. In addition, the integration in 
the polar direction is restricted in the range of =< ^ < tt, 
where Oi is given by equation (23). 

The numerical result is shown in figure 9. In figure 9, 
3/, the three times of the Eddington factor f{l3^—d/3/dr)^ 
is plotted by a contour map as functions of the velocity Po 
(abscissa) and the velocity gradient —d(3/dr\o (ordinate). 
Reflecting the breakup condition in the linear regime (see 
figure 6), the value of 3/ quickly drops at the outer edge 



of the one-tau photo-oval. Within the one-tau photo-oval 
region, the relativistic variable Eddington factor is well 
fitted by equation (33), except for the boundary region. 

4. Discussion 

Moment equations for relativistic radiation transfer 
have been derived in several literature (Lindquist 1966 
Anderson, Spiegel 1972; Hsieh, Spiegel 1976; Thorne 1981 
Thorne et al. 1981; Udey, Israel 1982; Schweizer 1982 
Flammang 1982, 1984; Nobih et al. 1991, 1993; Park 
2001, 2006 for spherically symmetric problems; Takahashi 
2007 for the Kerr metric). A complete set of moment 
equations for relativistic flow is given by the projected, 
symmetric, trace- free formalism by Thorne (1981). Since 
a moment expansion gives an infinite set of equations, one 
must truncate the expansion at the finite order by adopt- 
ing a suitable closure assumption, in order to make the 
transfer problem tractable. When one truncate at the 
second order, for example, the Eddington approximation 
(2) has been usually adopted as a closure relation. 

Such radiation moment formalism is quite convenient, 
and it is a powerful tool for tackling problems of rela- 
tivistic radiation hydrodynamics. However, its validity is 
never known unless fully angle-dependent radiation trans- 
fer equation is solved. 

Actually, the pathological behavior in relativistic radi- 
ation moment equations has been pointed out and exam- 
ined (Turoha, Nobih 1988; Nobih et al. 1991; Turoha et al. 
1995; Dullemond 1999). Namely, the moment equations 
for radiation transfer in relativistically moving media can 
generally have singular (critical) points. As a result, so- 
lutions behave pathologically in regions of strong velocity 
gradients. The appearance of singularities is explained be- 
cause we approximate the full transfer equations with a fi- 
nite number of moments (Dullemond 1999). For example, 
in the one-dimensional relativistic radiation flow using the 
Eddington approximation (2), where the moment equa- 
tions are truncated at the second order, the singularity 
appears when the flow velocity becomes /3 = v/c = ±l/ \/3 . 
Hence, under the traditional Eddington approximation, 
we cannot obtain solutions accelerated beyond c/V^, al- 
tough there exists a decelerating solution (Fukue 2005). 

As already discussed in the introduction, the invalidity 
of the Eddington approximation in such a relativistic flow 
can be understood as follows. In adopting the Eddington 
approximation (2), we assume that within the mean free 
path of photons the radiation field is isotropic in the co- 
moving frame. However, in the relativistic regime, where 
the velocity gradient becomes large and there exist the 
Doppler and aberration effects of photons, the isotropy of 
the radiation field may break down even in the comoving 
frame. 

Roughly speaking, the photon mean-free path £ in the 
comoving frame is l/{f<ip)^ where is the opacity mea- 
sured in the comoving frame and p is the proper density. 
When there exists a velocity gradient, dv/dz, the velocity 
increase at the distance of £ is estimated as 
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dv 
npdz 



dv 
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(35) 



In order for the radiation fields to be isotropic in the co- 
moving frame, this velocity increase should be sufficiently 
smaller than the speed of light. Otherwise, the usual 
Eddington approximation in the comoving frame would 
be violated. In addition, photons suffer from the Doppler 
and aberration effects, which change the radiation fields 
in the comoving frame. In such a case, we should modify 
the closure relation in the case of subrelativistic to rela- 
tivistic regimes, as in the case of optically thick to thin 
regimes. 

In order to improve the situation we are confronted 
with, instead of the usual Eddington approximation, we 
can adopt a variable Eddington factor, which depends on 
the flow velocity /3 {=v/c) and the velocity gradient d/3 / dr 
as well as the optical depth r (Fukue 2006; Fukue, Akizuki 
2006a, b; Akizuki, Fukue 2007; Koizumi, Umemura 2007). 
That is, in one-dimensional flows the variable Eddington 
factor f{r,/3,d/3/dr) is generally defined as 



P,o = f{T,/3,dp/dT)E,, 



(36) 



where Eco and Pco are the radiation energy density and 
the radiation stress tensor in the comoving frame, respec- 
tively. The closure relation in the inertial frame for one- 
dimensional flows then becomes 



cP(l - = cE{f - /3^) + 2F/3{1 - /), 



(37) 



where E, E, and P are the radiation energy density, the 
radiative flux, and the radiation pressure in the inertial 
frame, respectively (Kato et al. 1998, 2008). 

The function f{r,P,dP/dr) must reduce to 1/3 in the 
non-relativistic limit of = or in the uniform flow of 
d(3/dr = 0, whereas it would approach unity in the ex- 
tremely relativist ic limit of /3 = 1. Furthermore, in the 
sufficiently optically thick regime this function approaches 
1/3, while in the optically thin limit it reduces to an ap- 
propriate form determined by the geometry under the con- 
siderations. In the plane-parallel case, for instance, the 
variable Eddinton factor in the optically thin limit is an- 
alytically derived as 



1 -3/5 -h 3/52 



.-3/5 + /?2 ' 

since cE = 2E = 3cP in equation (37) in the optically thin 
limit (see also Kato et al. 2008; Koizumi, Umemura 2007). 

In these recent studies (Fukue 2006; Fukue, Akizuki 
2006a, b; Akizuki, Fukue 2007; Koizumi, Umemura 2007), 
they considered mainly the dependence of the variable 
Eddington factor on the flow velocity. Under the discus- 
sion above, in this paper we consider the dependence of 
the variabel Eddington factor on the velocity gradient, 
using the approach of the comoving observer's view point. 

As a result, for example, we derived the form of the rela- 
tivistic variable Eddington factor (33) and (34), where the 
Eddington factor slightly decreases for the subrelativsitic 
flow with the velocity gradient. In the extremely linear 
regime, the comoving radiation fields including P^^ and 
pyy perpendicular to the flow, and P^^, are calculated as 



cEr., 



pxx ^ 
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47r/o 1+ --^ 



' 47r/oi 
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5 dr 
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Using these results, the Eddington tensors become 
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^ 1+ 



d/3 
15 dr 
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15 dr 
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15 dr 
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(38) 

(39) 
(40) 
(41) 

(42) 
(43) 
(44) 



which satisfy the relation /^^ + /^^ + /^^ = 1. 

Thus, the Eddington factor in the direction of the flow 
slightly decreases, whereas that in the direction perpen- 
dicular to the flow slightly increases due to the velocity 
gradient. Since we have assumed that the emitted inten- 
sity is isotropic in the sufficiently optically thick regime, 
there is no aberration effect, and this result is due to the 
redshift effect, where the one-tau photo-oval is seen to 
expand for the comoving observer due to the velocity dif- 
ference. 

5. Concluding Remarks 

In this paper, from the viewpoint of the comoving ob- 
server we have analytically examined the Eddington factor 
in the relativistic flow accelerating in the vertical direc- 
tion; we have introduced the one-tau photo-oval observed 
by the comoving observer, and then calculated the comov- 
ing radiation fields and the Eddington factor. We thus 
have shown that the Eddington factor depends both on 
the flow velocity and the velocity gradient in the relativis- 
tic regime. 

Under the present plane-parallel case, the velocity in- 
crease in the downstream direction causes the density de- 
crease in the downstream direction; this effect is essential 
for the shape of the one-tau photo-oval and vessel. In the 
spherical flow, on the other hand, there exists a geomet- 
rical dilution effect, and the situation can be somewhat 
different. Such a spherical case is another problem. 

In this paper, in order to analytically examine the rela- 
tivistic variable Eddington factor and explicitly demon- 
strate its dependence on the velocity gradient, we re- 
stricted ourselves in the linear approximation, where the 
velocity can be expanded, and the velocity gradient is suf- 
ficiently small. However, the linear expansion of the ve- 
locity may not be inadequate when the velocity becomes 
large to be on the order of unity. Hence, the present ap- 
proximation would violate in the region of P ^ 1. We 
may examine such an extremely relativistic case under 
the same approach. 
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Moreover, in order to treat the problem in the sim- 
plest form, we only consider the sufficiently optically thick 
regime, where the one-tau photo-oval always exists. In 
the optically thin regime or in the large velocity gradi- 
ent case, there may appear the photo-vessel, where the 
optical depth becomes less than unity in the downstream 
direction. We should examine these general cases in the 
future. 

In the problem of line-driven winds with large ve- 
locity gradients, the Sobolev approximation is often 
adopted (e.g., Peraiah 2002; Castor 2004). The Sobolev 
method was also extended to the relativistic regime (e.g., 
Hutsemekers, Surdej 1990; Hutsemekers, Surdej 1993; 
Jeffery 1993, 1995a, b). Such a relativistic Sobolev ap- 
proximation is also local similar to the present case. 
However, it is not so simple to compare between the 
present continuum transfer case and the Sobolev's line 
transfer one. Whether the present approach can be ap- 
plied to the line transfer case, or the comparison with the 
Sobolev method is also the future work. 

Although the present treatment is limited in the lin- 
ear regime, we first examined the relativistic variable 
Eddington factor in the relativistic radiative flow under 
the physical situations, and obtained its dependence on 
the velocity and the velocity gradient in the explicit form. 

The relativistic variable Eddington factor for relativis- 
tic radiative flows discussed in the present paper is fun- 
damentally important in various aspects of relativistic as- 
trophysics with radiation transfer; i.e., black- hole accre- 
tion flows with supercritical accretion rates, relativistic 
jets and winds driven by luminous central objects, rela- 
tivistic explosions including gamma-ray bursts, neutrino 
transfers in supernova explosions, and various events oc- 
cur ed in the proto universe. 

The author would like to thank M. Umemura, T. 
Koizumi, and C. Akizuki for enlightening and stimulated 
discussions. This work has been supported in part by 
a Grant-in- Aid for Scientific Research (18540240 J.F.) of 
the Ministry of Education, Culture, Sports, Science and 
Technology. 



Hutsemekers, D., &, Surdej, J. 1990, ApJ, 361, 367 
Icke, V. 1989, A&A, 216, 294 
Jeffery, D.J. 1993, ApJ, 415, 734 
Jeffery, D.J. 1995a, ApJ, 440, 810 
Jeffery, D.J. 1995b, A&Ap, 299, 770 

Kato, S., Fukue, J., & Mineshige, S. 1998, Black-Hole 

Accretion Disks (Kyoto: Kyoto University Press) 
Kato, S., Fukue, J., & Mineshige, S. 2008, Black-Hole 

Accretion Disks - Towards a New Paradigm (Kyoto: Kyoto 

University Press) 
Koizumi, T., & Umemura, M. 2007, submitted to MNRAS 
Lindquist, R.W. 1966, Ann. Phys., 37, 487 
Mihalas, D. 1970, Stellar Atmospheres (San Francisco: W.H. 

Freeman and Co.) 
Mihalas, D., & Mihalas, B.W. 1984, Foundations of Radiation 

Hydrodynamics (Oxford: Oxford University Press) 
Nobili, L., Turolla, R., & Zampieri, L. 1991, ApJ, 383, 250 
Nobili, L., Turolla, R., & Zampieri, L. 1993, ApJ, 404, 686 
Park, M.-G. 2001, JKAS, 34, 305 
Park, M.-G. 2006, MNRAS, 367, 1739 

Peraiah, A. 2002, An Introduction to Radiative Transfer: 
Methods and applications in astrophysics (Cambridge: 
Cambridge University Press) 

Rybicki, G.B., & Light man, A. P. 1979, Radiative Processes in 
Astrophysics (New York: John Wiley & Sons) 

Schweizer, M.A. 1982, ApJ, 258, 798 

Shu, F.H. 1991, The Physics of Astrophysics Vol. 1: Radiation 

(California: University Science Books) 
Takahashi, R. 2007, MNRAS, in press 
Thorne, K.S. 1981, MNRAS, 194, 439 

Thorne, K.S., Flammang, R.A., & Zytkow, A.N. 1981, 

MNRAS, 194, 475 
Turolla, R., Nobili, L. 1988, MNRAS, 235, 1273 
Turolla, R., Zampieri, L., & Nobili, L. 1995, MNRAS, 272, 625 
Udey, N., & Israel, W. 1982, MNRAS, 235, 1273 



References 



Akizuki, C, & Fukue, J. 2007, submitted to PASJ 

Anderson, J.L., Spiegel, E.A. 1972, ApJ, 171, 127 

Castor, J. 2004, Radiation Hydrodynamics (Cambridge: 

Cambridge University Press) 
Chandrasekhar, S. 1960, Radiative Transfer (New York: Dover 

Publishing, Inc.) 
Dullemond, CP. 1999, A&A, 343, 1030 
Flammang, R.A. 1982, MNRAS, 199, 833 
Flammang, R.A. 1984, MNRAS, 206, 589 
Fukue, J. 2005, PASJ, 57, 1023 
Fukue, J. 2006, PASJ, 58, 461 
Fukue, J., & Akizuki, C. 2006a, PASJ, 58, 1073 
Fukue, J., & Akizuki, C. 2006b, PASJ, 58, in press 
Fukue, J., Kato, S., & Matsumoto, R. 1985, PASJ, 37, 383 
Hsieh, S.-H., Spiegel, E.A. 1976, ApJ, 207, 244 
Hutsemekers, D. 1993, ApJ, 417, 97 



